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In this paper we study an antiferromagnetic Heisenberg 
S — 2 quantum spin chain at T = with both interaction 
and on-site anisotropy, 



While the phase diagram of this chain for 5=1 and the 
isotropic point for S — 2 are by now excellently understood, 
contradictory scenarios exist for the S = 2 anisotropic phase 
diagram, which imply completely different mechanisms of the 
emergence of the classical S — > oo limit, hence the importance 
of the question. One main scenario predicts the emergence of 
a cascade of phase transitions not seen in the S = 1 case, 
based on an analysis using the Affleck-Kennedy-Lieb-Tasaki 
model. Another scenario is in favor of an almost classical 
phase diagram for S = 2, from which the S = oo phase dia- 
gram evolves quite trivially; the S = 1 case then emerges as 
very special with its dominant quantum effects. Numerical 
studies have so far not been conclusive. 

Using the Density Matrix Renormalization Group, we 
study gaps and correlation functions in the anisotropic S — 2 
chain. The question of the existence of hidden topological 
order in the anisotropic S — 2 chain is emphasised, as it dis- 
tinguishes between the proposed scenarios. Careful extrapo- 
lation allows us to show that the topological order is zero in 
the thermodynamical limit in all disordered phases, in par- 
ticular in the new phase interposed between the Haldane and 
large- D phases. The decay of the hidden order is substantially 
faster than that of S + S~ correlations and excellently fitted 
by power-laws typical of massless phases. This excludes the 
AKLT-model based scenario in favor of an almost classical 
phase diagram for the S — 2 and thus higher spin chains. 
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I. INTRODUCTION 

Interest in quantum spin chains as examples of low- 
dimensional strongly correlated quantum systems was 
greatly revived in 1983 by Haldane's famous conjectureEJ 
of a fundamental difference between half-integer and in- 
teger quantum spin chains, a phenomenon without equiv- 
alent in higher dimensional systems. 

In the meantime, Haldane's conjecture has been firmly 
established, while not rigorously proven. The numerical 
values of the predicted gaps and finite correlation lengths 
for isotropic integer spin chaias are by now extremely well 
known both in the S = MU and the S = 20119 cases, 



which constitute the extreme quantum limit. Further- 
more, for both cases Haldane's predictions of a massijie.- 
relativistic dispersion relation have been confirmedcrtjIiJJ. 

Haldane's conjecture was derived from a non-linear 
quantum sigma model, which is an effective low-energy 
continuum theory of the isotropic Heisenberg model in 
the semiclassical limit S — » oo. It does not provide im- 
mediate insight in the ground state of integer spin chains. 

The Affleck-Kennedy-Lieb-Tasaki (AKLT) modelHl 
provides a profound physical understanding of gapped 
spin liquids. The AKLT model is the Hamiltonian which 
gives rise to the following exact ground state: each spin of 
length S is decomposed into 2S completely symmetrized 
S = i spins. Pairs of such spins on neighboring sites are 
linked up by singlet bonds, giving rise to S singlet bonds 
between sites. In the case of S = 1, all important fea- 
tures of the isotropic Heisenberg model and the AKLT 
model, which has the Hamiltonian 
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are in qualitative agreementc 

A particularly striking feature of both models is the 
existence of a hiddcn—rjar^cjical topological order, also 
known as string order 



0(i,j) = (S?e i 



(2) 



which does not vanish in the limit \i — j\ — > oo. For the 
AKLT model, the value is |, for the Heisenberg model 
0.37432oEJ. In the case of a S — 1 chain, it describes a 
hidden antiferromagnetic order of Sf = ±1 spin states, 
which is diluted by arbitrary numbers of interposed Sf = 
spin states and thus not seen in the antiferromagnetic 
order parameter. 

The existence of the hidden order can further be asso- 
ciated with tho-breaking of a Zi x Zi symmetry in the 
Haldane phasetj, as well as with the existence of a four- 
fold degenerate ground-state of an open Heisenberg chain 
in the Haldane phaseEl For £? = 1, the success of the 
AKLT-model can be understood as it is a quantum dis.- 
order point, a remnant of a classical phase transitions. 
For higher spins, as we will see, both the existence of 
hidden order and its signification are not clear. 

In this paper, we want to make progress towards under- 
standing how the classical phase diagram emerges out of 
the substantially differeaJ, S =j=pL. phase diagram. Exam- 
ining different scenariosLnl3£jLJ, we will see that their 
fundamental difference regards whether or not there ex- 
ists hidden order in the S = 2 phase diagram for suitable 
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anisotropics. Our results will allow to decide quite clearly 
in terms of a scenario without hidden order for S = 2. 
We finish by addressing briefly why the AKLT model is 
not successful in explaining the anisotropic S — 2 phase 
diagram. 



II. S = 1 AND S = oo PHASE DIAGRAMS 

In real physical systems, such as NENP as the 
quintessential S = 1 spin chainEil, anisotropics are al- 
ways present, in the case of NENP an easy-plane on-site 
anisotropy D « 0.18E 2 ]. In order to link experiments to 
theory, it is therefore necessary to understand the phase 
diagram with anisotropies. At the same time, this may 
provide deeper theoretical insight into the interplay of 
the different phases. 

S = 1 phase diagram. The S — 1 phase diagram, as 
the physically most relevant one with at the same time 
the strongest quantum effects for an integer spin chain, 
has been studied extensively euad is by now extremely 
well understoo jyoy-NidO'Ea. Let us summarize the 
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main features (Figure |l]) . 

Both for substantially strong J z anisotropy and 
sufficiently small D, one finds Ising-like ordered 
(anti) ferromagnetic phases, depending on the sign of J z . 
The transition lines are not very far from those in the 
classical limit (see below). 
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FIG. 1. T — phase diagram of the S — 1 anisotropic 
Heisenberg antiferromagnet. Transition lines are schematic; 
the XY-Haldane transition line is at J z = exactly; the 
FM-XY transition line includes D = 0, J z = -1. For J z = 1, 
the Haldane-large-D transition is at D = 0.99(1); for D = 0. 
the Haldane- AFM transition at J z » 1.2. The XY phase 
actually comprises two different XY phases 1 and 2. 

All other phases are disordered, but otherwise physi- 
cally quite distinct: For large D, all Sf ^ states are 
suppressed and a singlet phase emerges, with a gap pro- 
portional to D, as the ga]i-is essentially the energy to 
excite one spin to Sf = ±lE3. Below that, there are two 
phases. Their separation line is strictly J z = 0, as first 
shown by den Nijs and RommelseES, and later on con- 



firmed numerically with increasing precisiorJl , - @. For 
J z > 0, there is a gapped quantum liquid, the Haldane 
phase, including the isotropic Heisenberg model, with 
the properties described above. For J z < 0, there is 
an XY-phase with power-law decay in the {S^~S~) cor- 
relations. In fact, the XY phase decomposes into two 
phases with different-jCorrelation functions, linked by an 
Ising-like transitionEatj. Phase transitions out of FM or- 
der are first order, from AFM order to the large- 1? phase 
first order, from Haldane to XY, and from XY to large- 
ly Kosterlitz-Thouless like, and the Haldane to large-Z? 
transition has been likened to a preroughening transition 
in RSOS modelsB. 
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FIG. 2. T = phase diagram of the classical anisotropic 
Heisenberg antiferromagnet. All transition lines are exact. 

Classical phase diagram. On the other hand, the T = 
phase diagram in the semiclassical limit S — > oo is very 
easily obtained (Figure [2]) . The quantum spins of length 
S are replaced by unit vector spins which do commute. 
This is in order as all terms in the Hamiltonian are of the 
order S* 2 , such that this factor can be scaled out. One 
finds three phases: 

1. J z > and D < J z — 1: ordered antiferromagnetic 
Ising-phase (AFM) 

2. J z < and D < — J z — 1: ordered ferromagnetic 
Ising-phase (FM) 

3. D > \ J Z \ — 1: ordered antiferromagnetic XY phase 
(XY): all spins in-plane 



III. EMERGENCE OF THE CLASSICAL LIMIT 

How does this classical phase diagram emerge from the 
5 = 1 phase diagram? In recent years, various scenarios 
have been proposed. 
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FIG. 3. Phase diagram for S > 1 as predicted by bosoniza- 
tion. The figure is not to scale. Adapted from Ref. 20. 

1. Bosonisation0 predicts for all integer spins 5 essen- 
tially the same scenario as for S = 1 (Figure |^) . Bosoni- 
sation can be expected to predict the nature of the phases 
found and their topology (i.e. which phases are joined by 
a transition), but not to yield quantitative values for the 
phase boundaries. 

For 5 = 1, bosonisation indeed predicts all phases 
and the correct topology of the phase diagram; it merely 
misses that the large- £>, AFM and FM phases meet re- 
spectively for very large anisotropics, beyond the limit 
of applicability of bosonisation. For higher spins, it is 
predicted in particular that for an isotropic exchange in- 
teraction J z = 1 there is exactly one transition from 
the Haldane phase to -the large-!? phase for D — > oo. 
From other argumentsEj it can be concluded that for 
large 5 the large- D phase will pushed out to D — > oo, 
whereas the gap of the Haldane phase vanishes; in the 
5 — > oo limit the gapless XY-phase emerges (but for the 
isotropic point, any D > will then favour in-plane spins 

Numerical studies0t3u9 are by now in excellent agree- 
ment that already for the 5 = 2 case there are not only 
two, but three phases when taking D — > oo for positive 
D and J z = 1 . While the results of Ref. E3are not recon- 
cileable with bosonisation, those of RefsT^ and [30] could 
be interpreted as being in accordance with bosonization, 
if one allows for a major deformation of the bosoniza- 
tion phase diagram. In the scenario we are going to show 
numerically below, bosonisation will be found to predict 
the correct phases, but we will also see that the necessary 
deformation of the bosonisation phase diagram changes 
it almost beyond recognition. .— . 

2. A second scenario was proposed by OshikawaO in 
the spirit of the Afneck-Kennedy-Lieb-Tasaki (AKLT) 
model (Figure Q). 

OshikawaHJ has shown that the hidden topological or- 
der parameter (the string order parameter) disappears 
exactly in the AKLT model, describing the isotropic 
chain, for all even integer spins 5 = 2,4,..., while it 
yields a finite value for odd integer spins 5 — 1,3,.... 
On the other hand, Oshikawa argues that the effect (as- 



suming for the moment an isotropic exchange interac- 
tion, J z = 1) of an increase in the easy plane anisotropy 
D > is to successively suppress the large S z states: first 
±5, then ±(5 — 1), and so on. This yields in his analy- 
sis a succession of spin chains with effective spin lengths 
5, 5 — 1, . . . , 1,0. The enviseaged succession of spin mod- 
els can be schematically drawn in the AKLT picture as 
in Figure |j| 




FIG. 4. T — phase diagram for S > 1 anisotropic Heisen- 
berg antiferromagnets proposed by Oshikawauj. Between the 
Haldane and large- D phases, S — 1 intermediate phases ap- 
pear, characterised by an alternatively finite and zero string 
order. 

As the string order parameter should be alternatingly 
zero and non-zero in the even and odd integer spin chains, 
Oshikawa predicts a succession of 5 + 1 phases with 5 
phase transitions, beginning with the isotropic Haldane 
phase (spin length S) and ending with the large- D phase, 
assimilated to a spin length 0. The transitions are char- 
acterized by the (dis) appearance of a finite string order 
parameter. The classical limit would then emerge by a 
softening of the quantum phase transitions and a disap- 
pearance of possible gaps in the 5 + 1 phases. 

The AKLT scenario works excellently for the 5 = 1 
case; and it is extremely well established— that a sim- 
ilar scenario, also proposed by Oshikawalia, of multi- 
ple phase transitions works extremely well to describe 
successive phase transitions in an isotropic Heisenberg 
integer spin antiferromagnet with increasingly strong 
dimerization of the interaction bonds, io-excellent agree- 
ment with jfield theoretical predictions^ and numerical 
verificationc3. All this makes this a very strong scenario. 

The existence of an intermediate phase, as should exist 
for 5 = 2, was first established in Ref. R_and later con- 
firmed in Refs. and [30| Oshikawa et al.Eil, using Quan- 
tum Monte Carlo on systems up to L = 160, claim to see, 
at least for D w 1.2 and J z = 1, a non- vanishing string 
order parameter, which they extrapolate to a thermody- 
namic limit of 1 x 10~ 3 at that particular point. They in- 
terpret furthermore the behavior of spin-spin correlations 
as indicative of a phase transition for J z rj 0. Further- 
more, a sign-alternation exists in the string order at the 
isotropic point, but not at D — 1.2. These results would 
support Oshikawa's scenario of multiple transitions. 
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FIG. 5. Mechanism proposed by Oshikawa for the emer- 
gence of the large-D phase out of the Haldane phase for J z > 
with increasing D. Circles represent spins of length S (here 
2), dots completely symmetrized spins- i. Bonds are singlet 
bonds. Up and down triangles symbolize symmetrized spins 
which contribute zero magnetisation and are not bonded. 

3. A third scenario results from numerical work by 
Schollwock et afl (Figure ||) . It actually has a .precursor 
in the works of Khveshchenko and ChubukovE3. They 
argue that the large-D phase moves as S 2 up to D — > 

00 for S — > oo, while the Haldane phase recedes to the 
isotropic point. The XY phase already present in the S = 

1 phase diagram invades the emptied part of the phase 
diagram; in the XY phase the hidden order parameter 
is zero in the thermodynamic limit. The 5 = 1 case 
with its strict J z = transition line between XY and 
Haldane phase emerges thus as particular. Numerically 
it is found that the Haldane phase is reduced to a small 
stripe next to the antiferromagnetic phase and squeezed 
out for S — > oo, as the gap is suppressed exponentially as 
e _7r ' s ; the large- D phase moves up to D — » oo, and the 
classical phase diagram emerges in a particularly simple 
fashion. There is no clear information on the behavior of 
the Haldane phase at its two end points. It is not clear 
whether (i) the Haldane phase touches the large- D phase 
at its top right end point and (ii) whether XY and AFM 
phase meet at its left end. If both are true, then not only 
the same phases, but also the same phase boundaries as 
in the bosonization picture emerge; the phase diagram 
looks still quite different. 

Numerically, the existence of an additional phase be- 
tween the Haldane and the large^D phase has been 
demonstrated!!! using the DMRgSM For J z = 1, the 
transition points were placed at D C1 — 0.04(1) and 
D C2 = 3.0(1). The intermediate phase was identified as 
gapless, with exponentially decaying (S Z S Z ) correlations, 
power-law like decay in the (Sf Sj) correlations, and a 
vanishing string order parameter, in support of the above 
picture. However, available computational ressources did 
not allow a very precise determination of the quantities 
involved (lengths up to 200, keeping 200 block states in 
the DMRG). 
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FIG. 6. T = phase diagram for S > 2 anisotropic Heisen- 
berg antiferromagnets proposed by Schollwock et al. The Hal- 
dane phase is almost squeezed out, the large-D phase moves 
up to infinity for large S. The XY phases extend rapidly with 
S, S = 2 is almost classical already. The transition line to 
large-D is not to scale (shown too low); the Haldane phase is 
exaggerated. 

Recently, Nomura and Kitazawa, combining exact di* 
agonalization and field theory techniques, have shownE2l 
that there are, on the J z = 1 line, two phase transitions, 
one from the Haldane phase to an XY phase, one from 
an XY phase to a large-D phase. They localise their 
transitions at D Cl = 0.043 and D C2 = 2.39. While the 
first result is in good agreement with DMRG results, the 
second is not. As the gap opens exponentially weak, the 
DMRG, which determines the gap directly, has difficulty 
in extracting such a transition point; fits to an exponen- 
tial law are hindered by the presence of small additive 
errors which become important when the gap approaches 
ps 0.001. However, evaluating new, more precise DMRG 
data with this fit formula shows that D C2 was overesti- 
mated in Ref. ^; with higher precision, we can show a 
gap down to D pa 2.5. In essence, the findings do agree, 
and the identification of the phase is the same in both 
works, at least on the J z = 1 line. The numerical data in 
Ref. |3(] leaves however open the possibility that the phase 
with a non- vanishing string order parameter proposed by 
Oshikawa might be located closer to the point where the 
large-D and antiferromagnetic phase meet directly. We 
will address this question now. 



IV. STRING ORDER IN THE S 
DIAGRAM 



2 PHASE 



1. Before we discuss our results on gaps and correla- 
tion functions, a remark on the calculation of the string 
order parameter is in order, as the crucial information is 
whether it is finite or not in the thermodynamic limit. 
Conventionally, the string order parameter is defined by 
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0(i,j) = {S?e™^ + i Si S*). (3) 

This means that the exponential string includes one of 
the two spins correlated. In the 5 = 1 literature, numer- 
ically often two alternative quantities are calculated: an 
order parameter with a string including both correlated 
spins, 

6 2 (i,j) = (S?e i * 1 ELi S l s . )) (4) 

and an order parameter including none of the correlated 
spins, 

0°(i,j) = (S?e i ^=* + i S *S?}. (5) 

Obviously, this changes some of the signs of contributions 
to the expectation value, depending of the actual value of 
the correlated spins. In the 5 = 1 case, it is obvious that 
a sign change occurs exactly in all contributions that are 
not zero, where both Sfj ^ 0. Therefore 

-O 2 (i,j) = O(i,j) = -O (i,j) (6) 

strictly. In the 5 = 2 case, this is no longer true: S z = ±2 
and S z = ±1 lead to different sign behavior in the three 
definitions. It is however, if computational ressources 
are stretched to their limit, more efficient for memory 
consumption reasons to calculate either O or O 2 in the 
DMRG method because then the string can be decom- 
posed into two symmetric parts, if i and j are chosen 
symmetric about the chain center of an open chain. Only 
one needs to be stored then. 

We have checked in lower precision calculations (200 
states kept, length 200) that in the region of importance 
the following inequality is always observed: 

\0 2 {i,j)\ < \0(i,j)\ < \O (i,j)\. (7) 

The very small difference in these correlations decreases 
with increasing \i — j\ both absolutely and relatively. 
Therefore, in the following, we always refer to the cal- 
culation of O (i,j). It provides both an upper bound to 
0(i,j) and converges to the same thermodynamic limit 
for i — j — > oo. In fact, on all our plots the difference be- 
tween the various correlations would be invisible to the 
eye. 

2. An observation that can be made numerically is that 
the string order parameter at D = 0, J z = 1 decays with 
alternating sign, whereas in the intermediate phase it de- 
cays with no sign-alternation. This might be attributed 
to a substantial difference between the various 5 = 2 
phases, as it has been done in Ref. |23|. Alternatively, this 
might naively be attributed to the fact that D > sup- 
presses both large S z values and leads to a very fast decay 
of S Z S Z correlations: With J z = 1, spin-spin correlations 
are antiferromagnetic for arbitrary D. This introduces a 
fundamental sign-alternation in the string order. What 
is effectively observed, now depends on the string: For 



D w and J z = 1, the close neighbourhood of the Ising 
AFM phase indicates that S z = ±2 states dominate in 
the string, contributing factors 1. The string order pa- 
rameter therefore exhibits sign-alternation. When D is 
successively increased, S z = ±1 and then S z = states 
dominate, introducing a sign alternation in the string for 
the former and no alternation for the latter case, such 
that the total order parameter is non-alternating and 
then alternating again. 

Calculating chains up to length L = 200 with the 
DMRG, we find for D > 0.02 that the string order pa- 
rameter shows sign-alternation for small distances, but 
no sign-alternation for longer distances. As we find (see 
below) a rather strong dependence of the string order 
parameter on the total system size if the string length 
becomes of the order of the system size, we have local- 
ized this "crossover" of the sign-behavior for various total 
system lengths. For D > 0.03, the crossover point is suffi- 
ciently small compared to the total system length that it 
may be considered almost converged; the numerical data 
does not exclude the possibility that the sign crossover 
can already be observed for smaller D < 0.02. We find 
crossover points as in Table Q. 

The sign crossover can be attributed to rather strong 
short ranged Neel like order, which is still present for 
small D, but does not survive for larger distances. It is 
of no importance whatsoever for larger values of D. The 
change in behavior is gradual and not directly linked to 
the critical value of D = 0.043 

For the large- 1? phase, the second change in sign be- 
havior is not observed up to D — 3.5, which is already 
deep in the large-D phase (gap of the order 0.05). 

This establishes that there is no direct relationship be- 
tween the phase transitions and the sign behaviour of the 
string order; a non-alternating decay is no unique feature 
of the intermediate phase. 

3. Furthermore, we have compared the decay of vari- 
ous correlation functions in the XY phase. In all points, 
we find a very fast exponential decay of S Z S Z correla- 
tions, and power-law decays in the string order and the 
5 + 5~ correlations. The decay of the string order is much 
faster than that of the typical XY phase correlations, 
even though the difference decreases towards the upper 
right corner of the XY phase. We see no change in be- 
havior between J z < and J z > 0. To give an example, 
decays follow laws x~ a , with a — 0.13 for J z = —0.5, 
D = 1.0 and a = 0.18 for J z = 0.5, D = 1.0 for the 
S+S- correlations and a = 2.11 for J z = -0.5, D = 1.0 
and a = 1.47 for J z — 0.5, D — 1.0 for the string order. 
This is indicative of the XY nature of the intermediate 
phase. 

4. Before we discuss results on the string order in the 
intermediate phase, let us first try to locate it with more 
precision. We consider only the part of the phase diagram 
with J z > and D > (Figure ||) . 

The transition line from the AFM phase to the Haldane 
or large- D phases can be easily located with the behavior 
of the S Z S Z correlation function. Considering chains up 
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to length L = 300, this transition line was determined 
with a precision of ±0.05 in the critical D-value. 

The Hakiane phase is very narrow (cf. former results 
for J z — iItl3), approximately 0.1 wide in units of D. A 
more precise determination has not yet been carried out, 
only at J z — 1.0 (considered above) and J z = 3.3 (see 
below) . 
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FIG. 7. Location of extrapolated gaps A for J z = 3.3 vs. 
anisotropy D. For D = 2.6, the chain is in the AFM phase. 
The bottom line A = 0.001 indicates our estimate of the 
precision limit. For smaller gaps, it cannot be clearly decided 
whether they are finite. The gapped phase on the left is the 
Haldane phase, the gapped phase on the right the large-D 
phase. Note that the gap curve would be almost flat, if the 
same scale were used for D and A. 

The crucial question is the location of the transition 
from the intermediate to the large-D phase, because the 
gap in the large- D phase opens exponentially slow. Here, 
we have calculated the gap in high-precision calculations 
(300 states kept, length 300) for various values of D for 
J z = 1; 2.5; 3.3. The fact that the gap is very small 
(or vanishes) implies an extremely large or divergent £, 
which has the advantage that a simple L~ ^extrapolation 
of the gap should provide a rather precise estimate of the 
gap (cf. the detailed discussion of such extrapolations 
m Refs. @,|). DMRG errors cancel to some extent as 
both ground and first excited state have similar trunca- 
tion errors (approximately 10 -8 or better for L = 600 
and M = 400 states kept). Variation of the number M 
of kept states reveals that we overestimate gaps slightly. 
Extrapolation in M as well as gap calculations for points 
where the gap is known to disappear let us estimate that 
gap values A > 0.001 are really finite. Consider Figure 
[71 For J z = 3.3, the truncation error is comparatively 
large, « 2 x 10 -8 . The very small gaps observed leave 



two scenarios: (i) between the Haldane phase (D « 2.65) 
and the large-D phase (for D larger than D w 3.0) there 
is a small intermediate phase; (ii) Haldane and large-D 
phase meet directly, at D ks 2.8. For smaller J z , the two 
phases can be clearly separated, still there is some un- 
certainty as regards the transition to the large-D phase. 
This leads us to the diagram shown in Figure the up- 
per transition line might in reality be somewhat lower 
(by no more than 0.2 units in D close to the upper right 
end of the XY phase; the transition point for J z = 1 is 
essentially exact). In any case, if we overestimate the 
size of the intermediate phase, there is no finite string 
order in the true intermediate phase if we don't find it in 
our oversized one. So our argument is not weakened. In 
fact, most of our points are far from the transition line 
anyways. 

5. Within the intermediate phase thus identified we 
are now going to calculate the thermodynamic limit of 
the string order for a large number of points. 

To discuss the possibility of an Oshikawa-like phase in 
the upper right corner, we will systematically investigate 
points for large J z . 

Our procedure to calculate the thermodynamic limit 
of the string order parameter is as follows. Using the 
DMRG, we calculate for a given set of (J Z ,D) the string 
order parameter 



0(N,L,M) = (S z e 



(8) 



with N = \i — j\, i, j symmetric about the center of an 
open chain of length L studied using a reduced Hilbcrt 
space of M states per block. 




FIG. 8. Location of points which were calculated for chain 
lengths L = 600, keeping 400 states for the DMRG. The tran- 
sition lines were obtained numerically; the AFM transition 
line is determined up to D = D c ± 0.05 or better; the Hal- 
dane phase has a width of « 0.05 or smaller and is thus hardly 
distinguishable from the XY phase, also due to the small gap. 
The D transition line is possibly slightly too high, as DMRG 
has difficulty in seeing extremely small gaps. 

To obtain the thermodynamic limit N — > oo, we first 
have to eliminate the DMRG typical dependence on M 
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of an order parameter for given N, L. This extrapola- 
tion can be done by fitting the string order against the 
DMRG truncation error. We find that for M = 400 the 
errors in the extrapolated results arc of no relevance to 
the following extrapolations. 
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O(N) = C(J z ,D)N- a(r ' D \ 



(9) 
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FIG. 9. Double-logarithmic plot of the string order 
0(N, L) for various system lengths L = 120, 180, 240, 300, 
360, 420, 480, 540, 600 for J z = 2.5, D = 2.3, with M = 400 
DMRG block states kept. The effect of the finite chain length 
as well as the convergence to a power law over several orders 
of magnitude are clearly visible. 

Having obtained 0(N,L) = limM->oo 0(N, L, M), we 
now have to take the thermodynamic limit first in L and 
N. The numerical result is an interplay of both the fi- 
nite N and L] as can be seen from the example shown in 
Figure |9|, the influence of finite L is strong. In Ref. |2^, 
these two limits are essentially taken at the same time, 
with N — L/2; to take this into account, a conformal 
field theory formula for correlations in systems with peri- 
odic boundary condition was used. As the power-law ex- 
tracted seems in very good accordance with our findings, 
which consider however much larger systems (L — 600 vs. 
L = 160), we attribute the deviations found in that ref- 
erence to the evaluation of the Monte-Carlo data used. 
Keeping N fixed, we find that the convergence of the 
string order parameter to thermodynamic limit value for 
fixed N is extremely well fitted by a L~ x law, if N is 
up to w 0.6L; for larger N, effects from the open chain 
ends start to dominate. For N less than « 0AL, the 
extrapolation correction is virtually negligible. 

We thus get 0(N) = limz-^ 0(N, L). This value 
can now safely be extrapolated to N — > oo. We find 
numerically that 0{N) vs. N can for N not too small 
(such that lattice effects disappear) be extremely well 
fitted by power laws such as 



as expected in a critical phase, where a varies slowly. 
The exponent a is given in Table || for our representa- 
tive highest precision points; all points calculated at less 
computational expense reflect the given results. The law 
expected in a massless phase from conformal field theory 
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for conformal weight h = h. We observe therefore a decay 
of the string order to zero in the thermodynamic limit in 
the intermediate phase, following a typical decay law of 
a massless phase, which we identify as an XY phase. It 
can also be observed that the exponent a decreases with 
increasing J z and increases with increasing D. 



V. CONCLUSION 

In our opinion, our numerical findings together with 
the analysis of Nomura and KitazawaEH show that the 
emergence of the classical limit is governed by the third 
scenario of an almost classical phase diagram for 5 = 2 
(Figure ||) . 

Why the failure of the AKLT model to account for the 
phase diagram of the S = 2 (and thus S > 2) phase 
diagram under the introduction of anisotropics, while it 
is extremely successful both for the isotropic Heisenberg 
model and for higher-S 1 spin chains with dimerization? 
In our view the failure is largely and simply due to the 
fact that the whole construction of the AKLT model(s) 
for various spin lengths S resides entirely on the use of 
singlet spin-spin bonds. At the isotropic point, each 
of these singlet bonds is separated from the neglected 
triplets by an energy roughly of the order of the gap 
energy. Dimerization is rotationally invariant and ac- 
counted for by singlet redistribution between neighbor- 
ing bonds. Anisotropies however destroy rotational in- 
variance, splitting the triplet states, and lower the ener- 
getic distance between the singlet and the lowest triplet 
state (s). A description of the chain simply in terms of 
singlets then quickly becomes inadequate. 

Obviously the S = 1 phase diagram has a very special 
property in that its XY and Haldane phases are sepa- 
rated by J z = strictly, which is linked to the fact that 
the Spil chain can be mapped to a six vertex times Ising 
mode£3, a mapping which is not possible for higher spins. 
The S — 2 phase diagram reveals no such special prop- 
erty. It might therefore be conjectured that the S = 1 
case itself is particular among integer spins and S = 2 
the really generic case. 
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TABLE I. Position of crossover from alternating to 
non-alternating sign behavior of the string order parameter 
for chain length L, J z = 1, and on-site anisotropy D. 



D L=50 L=100 L=150 L=200 

O00 = = = = 
0.01 

0.02 63 59 

0.03 41 39 35 

0.04 29 29 27 

0.05 27 23 23 23 

0.06 21 19 19 19 

0.07 17 17 17 15 

0.08 15 15 15 13 

0.09 13 13 13 13 

0.10 11 11 11 11 



TABLE II. Exponent of string order power law decay for 
various values of J z and D, extracted from chains of length 
L = 600, keeping M = 400 DMRG block states. 



r 


D 


a 


1.0 


1.5 


1.24 


1.8 


1.5 


1.14 


2.5 


2.3 


0.99 


2.7 


2.4 


1.00 


2.7 


2.7 


1.04 


3.0 


2.7 


0.95 


3.3 


2.8 


0.86 


3.3 


3.0 


0.90 


4.2 


3.6 


0.95 
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